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Abstract. In composed quantum systems, the presence of local dissipativc channels 
causes loss of coherence and entanglement at a rate that grows with the temperature 
of the reservoirs. However, here we show that if temperature is artificially added to 
the system, entanglement decay can be significantly slowed down or even suppressed 
conditioned on suitable local monitoring of the reservoirs. We propose a scheme to 
implement the joint reservoir monitoring applicable in different experimental setups 
like trapped ions, circuit and cavity QED or quantum dots coupled to nanowircs and we 
analyze its general robustness against detection inefficiencies and non-zero temperature 
of the natural reservoir. 
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1. Introduction 

Decoherence is considered the villain in quantum information processing, usually 
dragging previously entangled systems into separable ones. In the last few decades, 
different approaches to overcome its effects have been explored including quantum 
error correction codes [H El El Hj, dynamical decoupling [5], or encoding the system 
in decoherence free subspaces [6j [7] . 

Some other strategies make use of the decoherence process itself to try to protect 
entangled states. Reservoir engineering techniques, for example, explore the production 
of tailored made reservoir dynamics that have as steady solution the state to be 
protected [H [9]. Another option is to use the possibility to access the information that 
leaks to the environment to simply post-select favorable events, or to actively affect the 
system, as in the case of feedback methods [TQl EEH EES [13], [TU [15]. A common feature in 
these procedures is the requirement of some sort of non-local interaction to counteract 
the deleterious effect of decoherence on entanglement. These collective interactions can 
be easily created when the systems to be entangled are close together but can become 
a challenge for systems that are far apart. 

In the case of distant subsystems, the use of local strategies seems to be the most 
reasonable option from a practical point of view, although it is clear that any tactics 
based only on local operations would be useless to create entanglement. However, if 
an initial entanglement is shared between the distant parties, one could hope that local 
actions could reduce or even stop entanglement decay. 

Recently, local and non-local monitoring of quantum systems and its description 
in terms of quantum trajectories have been investigated in the context of entanglement 
decay [Ml [13, [TEl [191 ED] and local entanglement protection [211 E2]. In particular, 
Vogelsberger and Spehner [20] observed that the average entanglement could be 
preserved for systems coupled to local infinite temperature reservoirs under an optimal 
(in the sense that it maximizes the average entanglement) choice of a quantum trajectory 
unravelling. 

In this paper we present a proposal for the implementation of the above idea 
to preserve entanglement. First we derive the quantum dynamics conditioned on 
the measurement outcomes to show that the state evolves stochastically but remains 
entangled. Moreover, at any given time, the original state can be perfectly recovered 
as long as the detection signal is recorded with unity efficiency. Second, we show how 
to locally induce an artificial environment that reproduces the dynamics of decay and 
excitation reservoirs in order to mimic the infinite temperature condition. Finally, we 
provide a way to jointly, but locally, monitor both the engineered and the natural 
reservoirs in terms of quantum jumps or homodyne-like measurements. In this way, we 
give a physical interpretation in terms of a measurement prescription for the unravelling 
leading to entanglement preservation. Note that the local character of the proposed 
scheme allows its application to an arbitrary number of qubits, which makes it useful 
for different sorts of quantum communication and computation protocols. In the end, 
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we also discuss a specific experimental implementation of the scheme and analyse its 
robustness against imperfections. 

2. Preserving entanglement through local reservoir measurements 

We consider a composite system consisting of qubits that are weakly coupled to their own 
local and independent dissipative reservoirs at zero temperature with decay rates 7_ jQ . 
We also consider independently and locally engineered reservoirs that incoherently pump 
each subsystem with rates 7+, a . A practical way to produce these artificial environments 
will be presented in section [3j Under the Born and Markov approximations [23] the 
time evolution of the system is then given by the Lindblad [21] master equation: 

P = J^Ti.aZVi.oJP, (!) 

i.a 

where V[o~ ita ]p = a ita pa} a — l/2(a} a a ita p+pal a o~ ita ) with i = +,-, and a indicating the 
subsystem on which the operators are acting. <7_ a is the deexcitation operator for qubit 
a and <7 +jQ , its Hermitian conjugate. Note that in the particular case of 7+ < 7_ these 
engineered reservoirs simulate local thermal reservoirs of average number of photons 
^ = ~ ^-ry -i an d particularly of infinite temperature for 7_ = 7+ y]. 

This decoherence dynamics can be unravelled in terms of stochastic quantum 
trajectories in infinitely many ways, each one corresponding to a different scheme to 
continuously monitor changes in the reservoirs that act on the qubits [251 ESI E3 EEJ EH] . 
Quantum jumps, for example, describe the discrete clicks on a photodetector while 
homodyne-like measurements are represented by continuous stochastic equations. In 
what follows we will show how the local monitoring of ([1]) in terms of these kind of 
trajectories can preserve the entanglement in the system. 



2.1. Dynamics in terms of quantum jumps 

The most common unravelling is given by choosing the jump operators as the ones 
appearing explicitly in (JT]), i.e. Ji >a = \/ ^i, a dtOi^ a . The jump dynamics combines 
the action of randomly occurring quantum jumps given by the operators Jj ja with 
periods of no-jump evolution when the system evolves according to the operator 
Mnj = 1 — 1/2 Y2i a Ji a Ji,a- After the detection of a jump the system evolves to the state 

PJi a (t+dt) = T p" P j ^' (t)) wnere the denominator gives the probability (pi ta ) of detecting 
a jump i in the subsystem a. If no jump is detected in this interval, and assuming perfect 
detectors, the system then evolves to p^jCt + dt) = m njp(*)a-/ nj rp^g con ^j nuous 

monitoring of the system then indicates the sequence of operators that have to be applied 
at each time step determining, in this way, the trajectory of the system conditioned on 

I The infinite temperature regime is to be understood as the limit where n — ¥ 00 and 7 — > such 
that the rates 7_ and 7+ are the same and remain finite. In this way the system has equal chance of 
decaying or being excited. 
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that specific measurement record. In order to recover the unmonitored evolution of the 
density matrix given by the master equation, one has to average out all the jump and 
no-jump possibilities for the qubits: p(t + dt) = p^j p^j{t + dt) + ^2i a Pi,aPj ia (t + dt). 

To analyse the state dynamics, first note that when 7_ a = 7 + Q , = j a , the no-jump 
operator is proportional to the identity, 

M NJ = 1 - y ^2 la(o--, a O-+,a + <J+, a Cr-, a ) = 1 ~ "f^" ^ 
a a 

In this particular case, Mnj commutes with the jump operators Jj )Q ,, what means that 
the moment in which jumps are applied do not influence the final state of the system, 
only the number and type of jumps. In other words, the un-normalized conditional state 
of the system at time t for a given trajectory will be given by the sequential application 
of jumps to the initial state: 

p c (t) = (f[ j)j p(o) (n j)j , (3) 

where, once again, each J s corresponds to one of the possible J ija described above, and 
iV is the number of jumps detected in the particular trajectory. 

Also note that any of these jumps immediately destroys all the entanglement 
originally present in the system, since it identifies the state of the subsystem whose 
reservoir "clicked'. For example, if an spontaneously emitted photon is detected in 
the reservoir acting on a particular qubit (let's call it A), then the jump J-,a projects 
A into its ground state |0). Therefore this choice of monitoring scheme only allows 
the preservation of entanglement for trajectories where no jumps are detected, which 
happen with an exponentially decreasing probability in time. Therefore, the averaged 
entanglement over all possible trajectories also decays, even though it may still last 
longer than the entanglement of the unmonitored evolution [T9]. 

However, as mentioned before, this choice of jumps is not unique. In fact, new 
jumps defined as = Y^iUkiJi, with U (left) unitary, leave the master equation 
unaltered [jj while generating a distinct set of trajectories. Observe that the no-jump 
trajectory remains the same, independent of the unitary transformation U, and therefore 
the conditioned state at time t can still be found by the sequential application of the 
detected jump operators. 

We will now show that a particular choice of jump operators can in fact preserve 
the entanglement in every trajectory. If instead of a± >a we choose as new jumps the 
linear combinations 



Hadt hf a dt 



J y,» = l \l -^-(0--,a -0+,a) = y —^-Oy, a , (46) 

§ This unitary is not the most general transformation on the jumps that leave the master equation 
invariant. The addition of a complex number to the jumps with a corresponding change in the effective 
Hamiltonian is also possible |29j but we omitted it here for simplicity. 
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the jump operations will now be local and unitary, hence they cannot change the 
entanglement shared by the qubits. Since the product of Pauli operators is again a 
Pauli operator or the identity, the conditioned state p c {t) of the system at time t in any 
given trajectory will be given by 

Pc(t) = <ypauii,a p(0) ap au ii <a , (5) 

a a 

where each o-p au ii,a = EL ^s,a is one of {1, a Xia , o y ^ a Zja } and is determined by 
the sequence of jumps recorded on the respective subsystem a. In each trajectory 
entanglement is preserved because the system randomly jumps from one entangled state 
to another. However, provided that the parts keep track of their sequence of monitored 
jumps, the initial state can be recovered by locally reversing the final Pauli operation 
conditioned on the measurement results. Note that the scheme is valid for any number 
of qubits since it only relies on the fact that the jumps are local unitary operations and 
the state does not change between jumps. 



2.2. Dynamics in terms of diffusive unravellings 

The same idea can be applied to continuous in time stochastic unravelling. For perfect 
detection efficiency, the evolution of the system is given by [27J, EHl EH] 

dpy = ^ laT>[0-i,a]pY dt + ^2 ^ ~ fa.*) W d Ci, a + k.c] , (6) 

i,a i,a 

where d^ a are complex Wiener increments satisfying 

^ a (t)^(t) = ^, (7) 
d£i, a (t)d£* ta , (t) = dt u ijt J ao > , (8) 

with u being a complex symmetric matrix subject to the condition 1 1 2 < 1 for the 
matrix two-norm. Note that the 5 aa > indicates that all measurements are local and no 
correlation between the noises in the different subsystems are allowed in our model. 
The subscript Y in p is to remind that its evolution is conditioned on the measurement 
record given by the complex currents 

Y ija dt = ^ dtVl^i^a + Uij,aO-la) + dft,a- (9) 
3 

This time, the freedom in the choice of unravelling lies in the values of the noise 
correlations encoded in the matrix u. Each different u is associated with generalised 
homodyne-like measurements and corresponds to a particular way of mixing the signals 
from the monitored system and the noise introduced by external classical sources. To 
find the value that maximizes the average entanglement one can derive directly an 
equation for the entanglement in the system [30] as has been done for the case of zero 
temperature reservoir [TBI 120]. In what follows we will show that the choice Uij sa = — 1 
for i ^ j and Ua^ a = lead to perfect entanglement protection for a diffusive monitoring. 
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To see that, first note that with this choice of u the complex noises can be written 

as 

<%-,« = -j= , (10) 

, t ~dW ha + idW 2 , a 

d£+, a = , (11) 

with the real noises dW obeying dW^ a {t)dWj a '{t) = dt5ij5 aa >. Using these relations and 
C-,a = (o~x,a — io~y ta )/2, a +)OC = (o- Xja + io-y ta )/2 , we can express (jSJ) in a more convenient 
form: 

dpy = ^2 ~2~^t {o- Xl0i p?0- x , a + Vy,aP-?Vy,a ~ 2 Py) 
a 

~ J y [(^l.afv.a + ^2,0^,0),^] • (12) 

a 

Now it is easy to show that this evolution can be mapped into the application 
of local stochastic Hamiltonians defined as H a = y/^{dW\ ja a y + dW2, a o~ x ) / v2- The 
conditional state of the system at t + dt can be written as 

Py(t + dt) = e~ iHa p 9 (t)e iHa 

1 - iH a - ^f] p(t) (l + iH a - S) , (13) 



2 t 

where in the second step we expanded the exponentials up to first order in dt. Replacing 
the definition of H a and using the correlation properties of the real noises dW one obtains 
back Q12p . The evolution of the system therefore is equivalent to the application of a 
sequence of local Hamiltonians to the system, hence, unitary and local. Once again, the 
state evolves stochastically from one entangled state to another, now in a diffusive way, 
and entanglement is preserved for each trajectory. Moreover, because each subsystem 
has the record of their current signals and the local unitaries are state independent, 
whenever each party decide to use the qubits for any application, they just need to 
locally apply a final unitary matrix to revert the system back to its initial state. 



3. Engineering and monitoring the environment 

We proceed to show the engineering of the artificial temperature reservoir in atomic-like 
systems. Assume each qubit corresponds to the two lower levels of a three-level atom, 
|0) = \g) and |1) = \e) and that level |e) spontaneously decays into level \g) at a rate 
7_, emitting a photon circularly polarized to the right (annihilation operator or). The 
detection of this photon corresponds to a J_ oc o ge jump and upon its detection the state 
of the atom is certainly level \g). Let us also assume that levels \g) and \i) have the same 
m z projection of angular momentum and that an external 7r-polarized classical field of 
intensity Q is turned on at the \g) — > \i) frequency. In this case, if Q T then level \i) 
can be adiabatically eliminated and the overall effect of turning on the external field is 
to incoherently pump level \g) back to level |e) at a rate given by 7+ = and with an 
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Figure 1. Level scheme for engineering the reservoir. The qubit is encoded in levels 
\g) and \e), which can spontaneously decay back to \g). A classical field couples \e) 
and \i), which decays back to |e). When the decay rate T is much larger than the 
other frequencies of the problem, level \i) can be eliminated and an effective artificial 
excitation reservoir for the qubit is created. 



emission of a photon circularly polarized to the left (annihilation operator a^). Now, 
the detection of this photon corresponds to a J + oc a eg jump and upon its detection 
the state of the atom is certainly level |e). Note that if the same scheme is applied 
on all subsystems and the emitted photons are ignored, the dynamics of the qubits is 
described exactly by Furthermore, if ^ = 7_, then the incoherent pump, together 
with the spontaneous decay, simulate infinite temperature for the photonic reservoir. 

Within this model, it is easy now to describe the setup to implement the 
entanglement preserving monitoring. For each subsystem the scheme shown in figure [2^ 
has to be applied. If both photons have the same energy [jj] and polarized beam splitters 
(PBS) are placed before the detectors the photons will be combined in linear orthogonal 
polarizations an = {clr + cll)/V^ and ay = i(ai — au)/\/2. Their detection after the 
PBS will then correspond to the optimal jump operators given in (14 al) and (j-Xbj) . 

For the diffusive monitoring the scheme is shown in figure [2b. After the PBS 
the photons are combined with local oscillators at beam splitters and a homodyne 
measurement is performed. The currents obtained by subtracting the signals at the 
pairs of detectors D\ and and D3 and -D4 are 

h-idt = \J^{(aR ~ 4) - («£ - a R)) dt + dw ^ ( 14 «) 

I 3 . 4 dt = -iyJl((a R - 4) + (a L - a^ R )}dt + dW 2) (14&) 

with dWi real independent increments associated with the noise in the detectors. Now 
using the fact that the detection of photons a R and a^ correspond, respectively, to 
the action of cr_ and er + in the atomic system, one can combine the above signals 
to obtain the complex currents for the entanglement protecting unravelling, namely, 

j| in case uji e — uj eg = A then external fields can be applied to shift the electronic levels to produce 
photons of the same energy. 
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F_ = Ji_ 2 + H3-4 and Y + = — J 1 _ 2 + 2/3-4. Note that due to the symmetry of the 
reservoir, when one makes the correspondences — > <r_ and — > a + the average 
values terms in the signals vanish, indicating that for the protecting monitoring the 
signals contain only the noise parts. 

From a practical point of view, the difficult part in the above scheme can be the 
detection of spontaneously emitted photons due to the broad solid angle of emission. 
For this reason the scheme would be more suited to be implemented in one-dimensional 
systems where the emitters are embedded in broadband waveguides that allow the 
propagation of light in only one direction. Such systems are realizable nowadays in 
different experimental setups such as nanowires [3U [32] and superconducting circuit 
QED [331 El]- 



a) 



a R + a L 
-\wwv> • 




Figure 2. Monitoring scheme: photons corresponding to decay and excitation 
processes are collected and pass through a polarising beam splitter (PBS). For the 
quantum jump case (a), the detection of photons after the PBS will project the 
atomic state according to the optimal jump unravelling. For the diffusive case (b), 
the photons emerging the PBS arc combined with local oscillators and the signals 
recorded providing the currents in (|14q[) and (|146() . 

4. Analysis of imperfections 

In all the schemes here presented, the biggest problem is the detection inefficiency, 
specially because the added reservoir increases the emission rates in each qubit, therefore 
increasing the effect of the decoherence mechanism. We now proceed to calculate the 
effect of inefficiencies in our proposal. Inefficiency can be treated by including a rj 
parameter in the master equation, such that 77 = 1 means perfect detectors, and by 
rewriting (CQ) as 

P = V X] Ihc^WiAP + (1 _ v) Iho^^iAP- ( 15 ) 

The part proportional to 77 corresponds to the one about which we can obtain 
information, therefore the one we can unravel. The remaining part corresponds to a 
decoherence evolution, "slowed' by a factor (1 — 77), due to the information that is 
flowing to the reservoir and not being retrieved. 
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We shall now analyse the effect of detection inefficiencies in the quantum jump 
description. The most favourable unravelling described in Section 12.11 has a no-jump 
proportional to identity and jumps proportional do a x and a y . The difference now is 
that one has to add to the no-jump evolution a term corresponding to the decoherence 
part proportional to (1 — rj). Despite this difference, the analysis remains basically 
unaltered since the new no-jump evolution still commutes with the jump operations. 
It suffices to check this for one of the qubits since operations in different subsystems 
trivially commute. If one writes the local decoherence part in terms of o~ x and o y 
as in (|T2l one can easily verify that the action of a jump followed by decoherence 
(V[a x ](aipai) + V[a y ](aipai)) coincides with the evolution with decoherence followed by 
the action of a jump (ai(V[a x ]p + V[a y ]p)ai), where i is either x or y. 

As a result, the evolution of the density matrix of the system given that N jumps 
were detected will be given by the mixed state version of ©: 

Pc{t) = (^) o- PauU>a p(t) (g) a Pauliia , (16) 

a a 

where p(t) is given by the integration of equation p = (1 — rj) J2i a li,a < D[vi,<x[P- Note 
that this equation is equivalent to the unmonitored one except for the final application 
of local unitaries (that can be reversed since the parties know all the detected jumps), 
and the factor (1 — rj) corresponding to the detection efficiency. As previously noted, 
if the detectors are perfect, entanglement if fully preserved. Otherwise, the detection 
inefficiency becomes simply a multiplier factor on the decoherence rate. 

We are now in position to compare our scheme with both the monitored dynamics 
without the engineered reservoir and the unmonitored evolution. For that we consider 
the particular case of 2 qubits, A and B, that share initially a maximally entangled 
state of the form |^o) = ( 1 1 } + |10))/\/2. We also assume that the decay rates are the 
same (7,4 = 7b = 7). For this state, the unconditional entanglement, as measured by 
concurrence [35], decays as c(t) = e -7 * and c(t) = e~ 27< + e~ Alt /2 — 1/2 for zero and 
infinite temperature reservoirs, respectively [35] • The corresponding curves are shown 
by the solid lines in figure |3j black (a) for and red (b) for T . At zero temperature 
any local monitoring with 77 = 1 gives an average entanglement of e -7 * [TTJ HB] , coinciding 
with the master equation result (b). The concurrence for the entanglement conditioned 
on jump monitoring including the engineered reservoir can be calculated using (fTBl) 
and gives c(t) = e^ 1 "^ + e - 4 7(i-'?)*/2 - 1/2. The dashed blue curves (c) and (d) 
represent the cases with rj = 0.8 and rj = 0.9, respectively, and show that, for the time 
scale shown in the figure, monitoring preserves entanglement better than not adding the 
artificial reservoir. For any rj < 1, the dynamics will always have an infinite temperature 
component (with rate 7(1 — rj)) and the solution [36] will lead to disentanglement of 
the system at a finite time [37], which does not occur for initial Bell states at zero 
temperature. This means that the curves for the monitored system with artificial 
reservoir will have to cross the curve of the unmonitored system for the original zero 
temperature bath. The higher the detection efficiency, the longer our scheme will be 
advantageous as compared to the unconditional evolution. Obviously, for rj = 1 (curve 
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e) one has perfect protection for all times. Although this analysis has been presented for 
the case of two qubits, we should emphasise that it is also valid for multipartite systems, 
with perfect protection occurring for ideal monitoring and decline in performance as the 
detection inefficiency increases. 

For a diffusive unravelling, the inefficient evolution is not simple to calculate as 
in the jump case because now the sequence of evolutions do not commute with the 
unmonitored evolution. However we can still numerically simulate the equations and 
the result obtained follows exactly the curves for the jump monitoring. Therefore, since 
both kinds of monitoring perform equally well theoretically, choosing one monitoring 
scheme over the other is a matter of finding the one that can be implemented with the 
best possible detection efficiency 




Figure 3. Concurrence as a function of time for unmonitored reservoirs at infinite (a) 
and zero (b) temperatures. Dashed curves represent the average concurrence for the 
engineered reservoir under monitoring for r\ = 0.8 (c), 0.9 (d) and 1.0 (e). 

Finally, if the reservoirs have natural temperature, then the incoherent pump needs 
to be weaker in order to balance out both reservoirs. The equation in this case will be 
similar to the one with inefficiencies, with the natural excitation reservoir, which is not 
monitored, playing a role similar to that of an additional inefficiency to the detection 
scheme. 

In conclusion, we showed a way to engineer an artificial reservoir that increases the 
temperature of the system and yet can be used to decrease the rate of entanglement 
decay, or even stop it completely, when the system is properly (locally) monitored. With 
the use of simple optical elements, we showed how the monitoring scheme that preserves 
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entanglement can be implemented and derived a way to recover the initial state from 
the measurement record using only local unitary operations. In this way, the loss of 
information on the system due to the interaction with the environment, which, when 
irreversible, is ultimately the poison responsible for the decoherence process itself, can 
also be the vaccine used to fight decoherence as long as the information can be retrieved 
by a suitable monitoring scheme. 
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